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I.  INTRODUCTION 


Research  under  NASA  Grant  NsG-99-60  was  undertaken  as  a  part  of 
the  Jet  Propulsion  Research  Program  of  the  Department  of  Aeronautical 
Engineering  at  Princeton  to  "conduct  an  experimental-theoretical  study  of 
nonlinear  effects  in  combustion  instability  in  liquid  rocket  motors. 

This  basic  study  would  also  investigate  the  causes  and  magnitudes  of 
the  nonlinear  perturbations,  which  are  of  special  concern  in  the  design 
of  high  thrust  rocket  systems." 

II.  SUMMARY 

An  optical  model  of  a  no-flow  rocket  motor  has  been  completely 
designed  and  is  currently  in  an  early  phase  of  testing.  The  object  of 
this  optical  rocket  motor  is  threefold:  ( I )  to  establish  the  characteristics 
of  the  incoming  perturbation  -  in  the  current  system  a  perturbation  pro¬ 
duced  by  a  pulse  generator  utilizing  blank  cartridges  and  precision  diaphragms, 
(2)  to  aid  in  evaluating  the  velocity  probe  prior  to  its  insertion  into  a 
"hot"  rocket  motor,  and  (3)  to  observe,  for  the  applied  perturbation,  the 
shock  wave  patterns  produced  and  the  time  for  decay  under  varying  rocket 
motor  conditions. 

Preliminary  results  have  already  yielded  important  data  on  items  (I) 
and  (3)  and  under  item  (2)  the  probe  is  being  modified  for  considerably 
smaller  target  areas  because  of  the  shock  wave  strengths  encountered.  The 
shape  of  the  wave  as  well  as  its  strength  will  improve  the  basic  under¬ 
standing  of  the  problem  and  hence  influence  the  inputs  into  the  theoretical 
model  calculations. 

Nonlinear  rocket  tests  so  far  have  consisted  of  rechecking  one 
chamber-injector  combination  to  determine  the  stability  limits  prior  to  the 
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introductiofi  of  nonlinear  perturbations.  Actual  perturbation  tests  will 
begin  within  the  next  few  weeks. 

A  solution  to  the  theoretical  model  based  on  strong  shock  waves 
rotating  In  the  chamber  is  being  attempted  using  IBM  650  equipment.  A 
number  of  improvements  and  simplifications  have  been  made  In  the  progrem 
but  as  yet  no  solution  has  been  found. 

III.  DISCUSSION 

In  this  first  year  of  operation, the  effort  in  the  "Nonlinear 
Aspects  of  Combustion  Instability  in  Liquid  Propellant  Rocket  Motors" 
research  sponsored  by  NASA  has  concentrated  on  three  areas:  (I)  The 
optical  model  of  the  combustion  chamber,  (2)  the  design  of  hardware  and 
initiation  of  tests  on  the  nonlinear  rocket  motor  and  (3)  the  solution  of 
the  first  theoretical  model  relating  to  nonlinear  considerations.  The 
status  of  these  three  items  will  be  covered  in  this  report. 

A.  Optical  Model  of  the  Rocket  Motor 

The  optical  model  of  the  combustion  chamber  has  been  evolved  with 
a  number  of  objectives  in  mind.  In  order  to  study  the  effects  of  non- 
linearities  in  any  systerrv  one  must  first  be  able  to  measure  the  perturbation 
introduced.  In  previous  approaches  to  the  problem  only  empirical  relation¬ 
ships  were  sought.  The  onset  of  nonlinear  instability  was  measured  in 
terms  of  number  of  powder  grains  in  an  explosive  charge  (Ref.  I)  or  number 
of  psi  in  an  entering  stream  of  high  velocity  gas  (Ref.  2).  The  first 
portion  of  this  research  therefore  is  concerned  with  evaluating  exactly 
what  takes  place  at  the  station  where  the  propellants  are  injected,  since 
it  is  there  that  early  events  in  the  combustion  process  take  place. 

The  experimental  approach  to  evaluating  such  effects  within  the 
chamber  is  based  first  on  simulating  the  geometry  of  the  rocket  chamber  cross 


section  funder  no  through  flow  conditions]*  Using  optical  technigugs^, 
velocity  probes,  and  transient  pressure  Instriimentation,  the  time  history 
of  the  transverse  flow  conditions  due  to  the  entering  shock  pulse  are 
being  recorded.  Figures  I  and  2  Illustrate  this  experimental  apparatus 
located  in  the  newly  completed  laboratory  room  addition  to  the  nonlinear 
combustion  instability  rocket  cell. 

The  object  of  this  optical  rocket  model  is  threefold;  ( I )  to 
establish  the  character i sties  of  the  incoming  perturbation  -  in  the 
current  system  a  perturbation  produced  by  a  pulse  generator  utilizing 
blank  cartridges  and  precision  diaphragms,  (2)  to  aid  In  evaluating  the 
velocity  probe  prior  to  its  insertion  into  a  "hot"  rocket  motor,  and 
(3)  to  observe,  for  the  applied  perturbation,  the  shock  wave  patterns 
produced  and  the  time  for  decay  under  varying  rocket  motor  conditions. 

Concerning  the  experimental  apparatus  itself,  the  schematic 
of  Figure  2  is  most  helpful  in  the  overall  operation  of  the  system. 
Basically,  a  two-pass  schlleren  or  shadowgraph  system  Is  used.  A  spark 
is  fired  shortly  after  a  shock  pulse  enters  the  chamber  with  a  photo¬ 
graph  taken  of  the  resulting  density  gradients.  Details  of  the  individual 
items  of  apparatus  are  as  follows. 

The  optical  rocket  motor  is  identical  to  the  operating  rocket 
hardware  as  far  as  the  cylindrical  portion  of  the  combustion  chamber  is 
concerned.  In  both  cases,  the  cylindrical  section  is  four  inches  long 
and  is  composed  of  two  individual  chamber  sections,  one  containing  the 
pulse  gun  locations  (See  Figure  3).  Each  chamber  section  provides  for 
three  pressure  transducers  located  at  90^  intervals.  The  sections  may  be 
turned  with  respect  to  each  other  to  provide  pressure  instrumentation 
at  intervals  less  than  90°.  Thiswill  be  important  in  the  identification 
of  the  higher  modes  of  oscillation.  Although  It  Is  felt  that  the 


tangential  pulse  gun  location  (Figure  3)  will  receive  the  bulk  of  ex= 
perimental  attention,  a  radial  pulse  location  has  also  been  supplied. 

The  Injection  end  plate  in  the  optical  rocket  motor  is  designed 
to  provide  for  the  velocity  probes  at  the  same  locations  as  in  the  "hot” 
Injector  (Figure  4).  Ten  locations  are  provided.  In  the  case  of  the 
optical  injection  end  plate,  the  major  consideration  Is  the  optical  quality 
of  the  plate.  Manufactured  from  a  stainless  steel  plate  with  velocity 
probe  plugs,  the  face  was  ground  to  an  accuracy  of  +.  2  wave  lengths  of 
sodium  light.  Testing  to  date  has  caused  only  barely  detectable  errosion 
on  the  injector  plate  from  flying  bits  of  diaphragm  from  the  shots  fired. 
Tarnishing  and  smoke  deposits  are  removed  after  each  tests  using  alcohol 
applied  with  cotton  wadding.  In  the  absence  of  a  shock  pulse  small  sur¬ 
face  irregularities  from  the  grinding  process  are  visible  on  the  Injector 
face.  However  as  can  be  seen  in  figure  5  (which  describes  an  early  test 
on  the  optical  model  and  will  be  described  later  in  more  detail),  the 
shock  pulse  and  associated  turbulence  make  such  irregularities  unimpor¬ 
tant  in  the  interpreiation  of  the  optical  data.  Therefore,  grinding  to 
1/4"  wave  length  or  adding  a  nickel  coating  to  improve  reflectivity  and 
protect  the  surface  do  not  appear  necessary  at  this  time. 

At  the  other  end  of  the  chamber,  rather  than  the  convergent 
portion  of  the  nozzle,  there  is  instead  an  optically  flat  glass  plate.  In 
order  to  meet  the  requirement  of  chamber  pressures  as  high  as  1,000  psi 
with  a  nine  inch  chamber  diameter,  the  design  is  based  on  1-1/4  inch  thick 
glass  (the  thickest  in  normal  production).  For  the  higher  pressures 
tempered  glass  will  be  used  (tempered  glass  has  poorer  optical  qualities 
but  four  times  higher  strength),  while  ai  fiiu  lower  pressures  tf.e  polished 

plate  glass  with  an  ultimate  strength  which  allows  a  tOO  psi  (max.)  cham¬ 
ber  pressure  is  desirBtl<^.  Tiic  notiiod  Oi  clan  ping  should  provide  an 
additional  safety  factor  of  I.C.  AM  testing  is  tone  remotely,  in  any 


case,  due  to  tlie  uncertainness  In  the  glass  strengths. 

The  pulse  gun  Is  a  rnodifled  version  of  the  t/pe  used  by  /'iSrojet 
In  combustion  instability  tests.  The  principle  is  based  on  the  use  of 
a  carefully  sized  charge  of  rapidly  burning  gunpowder  to  burst  a  precision 
manufactured  diaphragm.  Thus  a  controlled  perturbation  is  supplied  to 
the  rocket  motor  being  tested.  The  pulse  gun  Is  shown  In  Figure  3  and 
consists  of  a  modified  Mauser  bolt  action  assembly  to  position  and  fire 
the  blank  cartridge.  The  solenoid  trips  the  conventional  trigger  mech¬ 
anism  which  in  turn  signals  the  spark  system  via  the  microswitch  located 
on  the  rear  of  the  bolt  action.  Internally,  as  is  shown  in  Figure  6,  the 
firing  pin  activates  the  primer,  the  flash  of  which  ignites  the  charge  of 
pistol  powder  causing  the  buildup  of  pressure  behind  the  diaphragm  which 
bursts  at  a  controlled  pressure  level.  Within  the  cartridge,  a  metal 
spacer  has  been  added  to  keep  the  dead  space  between  cartridge  and  dia¬ 
phragm  constant.  Three  grain  charges  are  currently  anticipated  -  45,  30 
and  15  grains.  In  the  present  tests,  only  the  15  grain  charges  have  been 
used  with  the  7500  psi  diaphragm. 

Returning  to  the  schematic  of  Figure  2,  the  next  item  is  the 
10  inch,  80  inch  focal  length,  optically  ground  mirror  1/4  wave  length). 

After  the  mirror  along  the  optical  path,  is  the  initial  light 
source,  the  knife  edge  and  the  photographic  equipment.  Light  i’s  supplied 
from  a  Spark  which  discharges  in  a  period  of  I  millionth  of  a  second 
(sufficiently  fast  to  reduce  possible  blurring  to  an  acceptable  level) 
with  energy  supplied  from  a  discharging  I 5Ky, .  5a/  condenser.  The  spark 
source  (borrowed  from  the  Gas  Dynamics  Group)  is  quite  simple  and  depend¬ 
able.  Initially,  the  optical  system  with  the  knife  edge  is  lined  up  using 
a  steady  source  of  light  from  a  zirconium  arc  lamp.  The  timing  of  the 
spark  is  of  particular  importance  since  the  shock  waves  degenerate  in  the 


order  of  a  few  milliseconds.  Therefore,  sn  electronic  timing  circuit 
activated  by  the  previously  mentioned  micro  switch  provides  the  appro¬ 
priate  time  prior  to  spark  discharge.  A  multiple-spark  firing  system, 
where  a  series  of  charged  condensers  are  discharged  in  turn  across  the 
spark  gap  will  be  utilized  shortly  in  the  test  apparatus.  Based  on  an 
automotive  distributor  principle,  where  rotational  speed  will  determine 
time  intervals,  problems  center  about  maintaining  a  small  tickle  charge 
so  that  at  either  the  distributor  points  or  the  spark  gap  no  problem  of 
a  nonlonized  condition  exists  and  thus  restricts  discharge. 

The  flash  of  the  pulse  gun  itself  has  been  shown  to  be  impor¬ 
tant  as  a  source  of  fight,  with  intensities  of  the  same  order  as  that  of 
the  spark.  Since  the  shutter  must  be  opened  in  the  darkened  test  room 
awaiting  the  spark,  filters  limiting  light  to  the  ultra  violet  range 
were  used  in  addition  to  Increasing  the  spark  intensity  to  reduce  flash 
1 nterf erence. 

Currently,  the  photographic  equipment  is  quite  simple  and  con¬ 
sists  of  a  remotely  operated  shutter  on  a  portrait  camera  where  the  image 
has  been  prefocused  on  ground  glass.  The  camera  is  off  to  side  due  to 
space  requirements  and  hence  a  flat  mirror  is  used  (See  Figure  I).  Since 
with  a  two-pass  schlieren  system,  alignment  must  include  an  angular 
deviation  in  the  overall  optical  path,  a  number  of  images  are  produced. 

At  the  knife  edge,  the  brightest  image  is  selected  and  the  remainder  are 
masked  out.  When  the  multiple  spark  system  is  used,  the  photographic 
complexities  will  increase  accordingly.  Since  requirements  are  such  that 
six  photographs  must  be  taken  in  an  approximately  three  millisecond  in¬ 
terval,  a  drum  camera  is  called  for.  A  6  inch  diameter  drum  camera  is 
avdilaole  and  will  require  rotational  speeds  of  12,000  RPM  to  separate 
tl]o  phoTojrapi  s  isifc  I tih.  BlJrring  during  the  period  of  the  spark 
would  bo  lirited  fo  approxlmofuly  .001  inch  and  v/ould  fe  quite  acceptable. 


The  Instrumentation  used  In  these  tests,  other  thar  the  optical 
system,  has  already  been  mentlonad  briefly,  namely,  the  pressure  trans¬ 
ducers  and  the  velocity  probes.  The  latter  require  additional  explanation, 
however,  since  they  must  still  be  considered  in  the  development  stage. 
Figure  7  gives  a  cross-sectional  representation  of  the  velocity  probe. 

An  extremely  simple  device.  It  Is  shown  to  consist  of  the  sensing  target, 
the  pressure  seal,  ball  linkage,  strain  gage  column,  second  pressure 
seal  and  eJectrical  connections.  The  strain  gages  are  supplied  with  either 
3  or  6  volt  DC  and  provide  a  I  millivolt  signal  (with  3  volts)  for  lO  pound 
load  on  the  target.  One  Important  area  of  development  has  been  in  the 
case  of  the  targets,  where  a  material  capable  of  withstanding  rocket 
temperature  (5000^F4),  end  high  heat  transfer  rates  (15  BTU/in^/sec  or 
more  are  present  with  strong  transverse  combustion  instability)  was  sought. 
Various  materials  were  tested  in  motors  operating  on  the  "linear"  program 
and  included  lucite,  graphite,  ceramic  coated  stainless  steel,  ceramic 
coated  copper,  boron  nitride,  astrolite  and  other  glass  reinforced  com¬ 
position,  and  GRB  silicon  carbide.  The  only  materials  which  withstood  the 
conditions  (from  minus  several  hundred  degrees  to  operation  for  two 
seconds  at  approximately  5000°F)  were  the  GRB  silicon  carbide  and  the 
ceramic  coated  copper  and  in  the  latter  case,  a  new  ceramic  coat  was 
necessary  for  each  test.  More  recent  tests  of  the  GRB  silicon  carbide 
target  are  shown  in  Figure  8,  which  illustrates  how  well  this  material 
stands  up  under  the  high  heat  flux  rates  as  well.  These  tests  were  run 
on  the  square  motor  in  the  "linear"  program  (Reference  3). 

Sizing  of  the  targets  is  under  further  consideration.  Testing 
In  the  optical  rocket  motor  has  shown  the  size  to  be  too  large  for  close 
proximi'^y  to  the  incoming  pulse  (the  darkened  spot  in  Figure  5  indicates 
this  position).  PI  Stic  end  aluminum  targets  were  damaged  in  the  early 
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tests  Indicating  forces  well  beyond  prellmlnar/  estimates.  More  recently  a 
hollow  stainless  steel  rod  was  used  with  no  damage  resulting.  Until  a  fuller 
knowledge  of  the  shock  wave  pattern  is  obtained,  a  probe  cannot  be  calibrated 
in  the  optica!  apparatus  directly.  However,  with  the  correct  target  shape 
known  calibration  will  take  place  in  available  shock  tubes. 

Since  the  optical  rocket  has  only  been  in  operation  over  the  past 
few  weeks,  any  data  must  be  considered  preliminary.  The  photograph  shown  in 
Figure  5  illi.-siratos  one  type  of  data  that  is  being  sought  with  this  apparatus. 
The  curved  shock  pattern  when  compared  to  the  point  where  the  shock  pulse  Is 
introduced  would  indicate  that  although  the  pulse  is  oriented  tangentially, 
the  shock  initially  travels  across  the  chamber.  The  point  of  injection  is 
seen  to  be  displaced  from  the  center  of  the  resultant  spherical,  wave. 

However,  port  shape,  wall  reflections  and  the  variation  In  gas  properties,, 
due  to  the  incoming  shock  pulse  could  account  for  this,  Sch Merer  photos 
at  various  points  during  the  shock  growth  should  help  answer  such  questions. 

Pressure  and  velocity  probe  data  from  these  tests  have  been  found 
to  add  additional  information  as  to  the  sequence  of  events  following  the 
first  pressure  rise  (which  is  difficult  to  observe  due  to  resonance  in 
the  transducers  in  the  presence  of  a  step  function  input).  The  pressure 
records  indicate  that  at  least  by  the  third  distinct  pulse,  pressure 
pulsations  are  traveling  tangentially  in  the  direction  of  the  pulse  gun 
orientation. 

B .  The  Nonlinear  Rocket  Motor 

Based  on  the  wide  variety  of  tests  performed  on  the  "linear" 


combustion  instability  program,  it  was  decided  to  use  an  injector-chamber 
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combination  that  had  already  exhibited  a  clear  set  of  combustion  stability 
limits  when  mixture  ratio  was  varied.  The  combination  chosen  was  the  500 
pound  thrust  level,  1.4  design  mixture  ratio  spuds  using  an  eight  Inch 
Injection  diameter  in  a  nine  Inch  diameter  chamber. 

The  planning  behind  the  testing  on  this  rocket  motor  Is  shown 
In  the  sctiematic  of  Figure  9.  Starting  from  a  known  set  of  stability 
limits  for  the  standard  rocket  conditions  (only  combustion  noise  present 
and  starting  transients  guarded  against  with  destructable  barriers. 

Ref.  4)  the  perturbation  level  will  be  raised  in  a  controlled  fashion  and 
any  shift  in  stability  limits  will  be  noted.  The  perturbation  will  be 
supplied  from  the  pulse  guns  tested  on  the  optical  rocket  motor.  It  is 
hoped  that  a  relationship  between  the  deviation  of  the  stability  limits 
can  be  correlated  with  the  basic  parameters  associated  with  the  rocket 
motor  design. 

Currently  the  standard  conditions  are  being  recfiecked,  since  the 
velocity  of  any  comparisons  rests  on  these  standard  limits.  This  same 
injector  type  has  also  been  involved  in  recent  tests  with  the  "linear" 
program.  Certain  starting  d  I  f  f  icul  ti  es,  which  have  resulted  In  the  loss-of 
rocket  hardware  (Ref.  3),  have  limited  the  number  of  tests  on  the  nonlinear 
hardware  until  the  di fficulty  could  be  cleared  up.  New  starting  procedures 
and  safety  instrumentation  are  now  incorporated  in  the  test  operation  of 
both  researches  and  testing  will  increase  rapidly  in  the  next  few  weeks. 
Since  this  test  stand  is  closely  tied  to  the  results  and  checkout  of 
instrumentation  on  the  optical  rocket  motor,  testing  on  both  stands  must 
proceed  in  a  coordinated  fashion. 

Ceriain  other  nonlinear  data  which  is  constantly  being  collected 
as  a  by-product  of  the  "linear"  program  will  also  be  checked  out  in  far 
greater  detail  on  fhe  nonlinear  rocket  motor. 


C.  The  Nonlinear  Theoretical  Wodel 


A  complete  description  of  the  nonlinear  theoretical  model  is 
Covered  In  Appendix  A  with  the  detailed  equations  stated  in  Appendix  B, 

From  the  theoretical  point  of  view,  a  knowledge  of  velocity  and 
pressure  perturbations  Is  Important  In  the  determination  of  inputs  into 
the  numerical  sol ut Ions  currently  in  progress.  Briefly,  the  theoretical 
attack  to  date  has  Involved  a  model  consisting  of  a  cylindrical  annular- 
type  chamber  with  a  constant  flow  of  premixed  gases  entering  at  one  end. 

The  resulting  theory  would  also  be  applicable  In  a  conventional  cylindrical 
chamber  near  the  wall. 

In  the  model  pressure  disturbances  rotate  in  the  transverse 
direction  with  constant  velocity  and  constant  wave  length.  .All  combustion 
Is  assumed  to  occur  when  the  premixed  gas  Is  swept  by  the  disturbance  the 
first  time,  thus  a  generalized  Chapman- Jouget  detonation  criterion 
applies.  The  release  of  energy  In  combustion  is  considered  the  driving 
mechanism  for  the  system. 

From  these  and  other  considerations  the  four  basic  algebraic 
equations  have  been  derived.  They  involve  (1)  the  Chapman-Jouget  condition, 
(2)  the  Hugoniot  relation,  (3)  the  I soenergetic,  and  (4)  isentropic 
relations  applied  between  the  disturbances.  Since  all  pressure  waves  ere 
assumed  to  be  purely  transverse,  the  axial  velocity  Is  constant  throughout 
the  flow  and  the  axial  momentum  equation  Is  trivial.  Important  factors 
in  these  conservation  equations  are  the  lengths  of  the  detonation  region 
and  the  shock  region  of  the  disturbance.  They  appear  as  weighting  factors 
in  the  equations. 

Acquisition,  (on  a  rental  basis)  of  an  IBM  1620  computer  late  this 
summer  will  greatly  alleviate  the  major  problem  In  the  solution  of  the 
theoretical  model  so  far,  1hat  of  having  adequate  machine  time.  Recent 


breakdowns  and  the  heavy  load  of  other  University  researcl-irava  limited 
access  to  the  550  machlra  to  only  a  fcv#  hours  per  week. 

Solution  of  the  first  theoretical  model  will  be  followed  by  the 
more  complex  case  where  the  curvature  of  the  shock  waves  is  -taken  Into 


Consideration 


/IPP£KDIX  A 

Theoretical  Model  of  Nonlinear  Combustion  Instability 

For  Initial  simplicity  in  determining  a  theoretical  model  for 
the  phenomenon  of  nonlinear  combustion  Instability,  a  cylindrical  annular- 
typo  Combustion  chamber  has  been  considered  with  a  constant  flow  of  pre¬ 
mised  gases 'entering  at  one  end.  The  resulting  theoretical  treatment 
may  also  be  applied  near  the  wall  of  a  conventional  cylindrical  chamber. 

In  the  model,  one  or  more  pressure  disturbances  (assumed  concen¬ 
trated  in  an  axial  planed  rotate  In  a  transverse  direction  with  constant 
velocity  and  constant  wave  length.  Each  disturbance  Is  Identical  In  nature 
and  extends  the  length  of  the  chamber.  All  the  combustion  occurs  when  the 
premixed  gas  Is  swept  by  the  disturbance  for  the  first  time.  In  the 
following  sweeps,  no  combustion  occurs,  just  non-lsentroplc  compression. 
Therefore,  the  portion  of  the  disturbance  near  the  Injector  Is  a  detona¬ 
tion  (assumed  to  be  the  Chapman-Jouget  type  since  this  Is  the  only  stable 
type  one  can  imagine}.  The  other  portion  of  the  disturbance  Is  a  pure 
shock  wave.  In  between  these  concentrated  compression  waves,  expansion 
waves  exist  and  these  must  move  in  a  transverse  direction,  also. 

The  release  of  energy  In  combustion  Is  considered  the  driving 
^  * 

mechanism  for  the  system. 

If  the  reference  frame  Is  fixed  to  this  disturbance,  the  problem 
becomes  a  two-dimensional  steady  type,  and  independent  variables  are  the 
axial  and  transverse  coordinates.  Radii  are  assumed  large  enough  to  neg¬ 
lect  centrifugal  effects  and  radial  variations  In  fluid  properties  are 
neglected.  In  the  formulation  of  the  theor>  average  values  are  used  for 
the  velocities,  pressures  and  densities  at  the  surface  of  the  concentrated 
disturbances.  Averages  over  each  of  three  regions  are  used.  These 


readers  are;  tfie  regfon  of  unburned  gases  at  th@  iipstreaui  surfae®  of  the 

detonation,  the  region  of  burned  gases  at  the  upstream  surface  of  the  shock, 

«  * 

and  the  region  of  burned  gases  at  the  downstraam  surface  of  the  concentrated 
disturbance. (See  Figure  |0), 

The  assumption  Is  made  that  the  average  pressui  es  are  the  same 
for  the  burned  and  unburned  regions  st  the  upstream  surface  of  the  dis¬ 
turbance,  All  coefficients  of  non-unI formity  which  would  appear  In  this 
formulation  are  assumed  to  be  unity. 

The  conservation  equations  for  mass,  energy,  and  the  transverse 
component  of  momentum  are  applied  across  the  surface  of  the  disturbance 
as  a  whole.  Since  all  pressure  waves  ere  purely  transverse,  the  axial 
velocity  is  constant  throughout  the  flow  and  the  axial  momentum  equation 
adds  no  additional  Information,  Important  parameters  appearing  In  these 
conservation  equations  are  the  lengths  of  the  detonation  region  and  the 
shock  region  of  the  disturbance.  They  appear  as  weighting  factors  In  the 
equations.  These  lengths  must  be  determined  as  functions  of  the  pressures, 
densities,  end  velocities  at  the  surfaces  of  the  disturbance  and  at  the 
Injector  surface.  This  is  accomplished  by  means  of  an  Investigation  of  ths 
particle  trajectory  in  the  unburned  gas  region. 

The  assumption  is  made  that  the  pressure  gradient  In  the  trans¬ 
verse  direction  Is  constant  near  the  Injector.  The  Isentroplc  condition 
Is  applied  throughout  the  unburned  gas  region,  and. allows. a  determination iof. 
the  functional  dependence  of  these  lengths  upon  certain  pressures,  den¬ 
sities,  and  velocitl  as. 

The  application  of  the  Isoenergetic  condition  in  this  region  of 
unburned  gasss  adds  a  new  unknown,  the  velocity  of  the  disturbance  relative 
to  the  chamber.  This  is  also  the  transverse  velocity  component  (in  the 
frame  of  reference)  of  the  gas  at  the  injector  surface.  The  isoenergetic 


^^3 

condition  may  be  eppiled  at  ths  end  of  the  analysis  to  determine  the 
velocity  of  the  disturbance. 

The  Isentropic  and  Isoenergetlc  conditions  are  applied  across 
the  region  of  expanding  burned  gases  to  relate  average  fluid  properties  at 
the  downstream  surface  of  one  disturbance  to  properties  at  the  shock  por¬ 
tion  of  the  upstream  surface  of  the  next  disturbance.  These  relations 
actually  connect  properties  on  one  side  of  the  disturbance  to  properties 
on  the  other  side  because  of  the  cyclic  nature  of  the  model. 

The  Chapman-Jouget  condition  is  applied  across  the  detonation. 

That  Is,  the  velocity  of  propagation  of  the  detonation  wave  Is  minimized 
with  respect  to  the  pressure  ratlo^across  the  disturbance.  Note  that  the 
thermodynamic  properties  In  front  of  the  detonation  are  considered  constant 
in  the  differentiation.  It  must  be  noticed,  however,  that  this  particular 
way  of  applying  the  Chapman-Jouget  condition  contains  a  certain  element  of 
arbitrariness. 

The  foil  owing  parameters  are  assumed  to  be  known  Initially; 
quantity  of  heat  released  In  combustion,  specific  heat  ratio  (assumed  the 
same  everywhere),  chamber  length,  wave  length  of  disturbance,  average 
thermodynamic  properties  of  premixed  gases  at  injection,  and  the  axial 
velocity  Component  of  premixed  gases  at  injection.  The  unknowns  are  the 
average  values  for  the  pressures,  densities,  and  transverse  velocities 
along  the  disturbance  and,  also,  the  propagation  velocity  of  the  disturbance 
relative  to  the  chamber.  There  are  now  equal  numbers  of  unknowns  and  equations. 

For  the  equations  themselves,  the  following  nomenclature  will  be 


0~  =  speed  of  sound 


\p  = 


pressure 

density 


used : 


£t= 

L- 

/- 

r  = 

\r 

Subscripts 

/tt- 

th  = 

f 

L  = 


vaJocity 
wav9  length 
chamber  length 

lerg+h  of  ffatonatlon  portion  of  disturbance 

length  of  shock  portion  of  disturbance 

ratio  of  specific  heats 

velocity  of  disturbance  relative  to  chamber 

longitudinal  velocity  of  gas  throughoilt  the  field 

energy  per  unit  mass  of  reactant  released  In  combustion 

conditions  at  unburned  region  at  upstream  surface  of 
d i sturbance 

conditions  at  burned  region  at  upstream  surface  of 
disturbance 

conditions  at  downstream  surface  of  disturbance 
conditions  at  surface  of  Injector 


It  is  convenient  to  nondimenslona I i ze  as  follows: 


Q  '=  fAf 


{= 
h  - 


The  isentropic  condition  is  applied  along  a  streamline  in  the 
unburned  gas  region  to  obtain  the  following  relation  between  conditions 
at  the  injector  and  at  the  surface  of  the  disturbance  u. • 
The  length  of  the  detonation  zone  ^  Is  easily  found  as  a  function  or 


0  and  2  by  relating  the  mass  flow  through  Ihe  Injector  to  the  mass  flow 


through  the  detonation  region.  Thus, 


so  that. 


The  three  conservation  equations  may  be  applied  across  the  dis¬ 
turbance.  Mass  and  momentum  equations  are  used  to  eliminate  Uj  and  from 
the  conservation  of  energy  equation.  This  results  in  an  equation  of  the  form 


The  Isoenergetic  condition  applied  in  the  region  of  burned  gases 
gives  another  relation  between  properties  at  the  surfaces  of  the  disturb¬ 
ances.  Again,  the  equations  of  the  conservation  of  mass  and  momentum  are 
used  to  eliminate  velocities.  The  relation  is  of  the  form 


The  isentropjc  condition  applied  in  the  region  of  burned  gases 
provides  still  another  relation. 


Q/f  =  (p 


This  equation  is  used  to  eliminate  from  Equations  (I)  and  (2) 


to  obtain 


C(‘P, 


The  Chaprran-Jouget  cooditJofi  Incficates  that 


Equations  (4),  (5)  and  (6)  are  three  simultaneous  equations  In 


the  variables  0  ,  0  ,  and  E  and  the  parameters  ^  »  A  ,  and  JT  . 
The  problem  Is  then  one  of  solving  these  three  equations  for  the  values  of 
0  ,  0  ,  and  E  for  giv'en  values  of  ^  #  A  and  X  .  The  three 

equations  are  given  in  full  form  in  Appendix  B, 

This  theoretical  problem  is  In  the  process  of  being  solved  using 
an  IBM  650  computer.  An  Iteration  scheme  Is  being  used,.  Rather  than 
arbitrarily  choosing  starting  values  for  0  ,  9  ,  and  Z  in  the  Iteration 
scheme,  theIBM65C  i s  ca Iculatl ng  the  values  of  the  functions  F(0,  9,  2) 
and  H(0,  9,  1)  at  each  of  a  lattice  of  points  in  0  ,  0  ,  E  space. 

From  this  plot  of  points  a  rough  estimate  of  the  geometry  of  the  F  =  0  , 

G  =  0  ,  and  H  =  0  surfaces  may  be  expected.  At  this  time,  the  values 
of  F,  G,  and  H  have  been  determined  at  a  lattice  of  points  covering 
the  space  I  0  25,  I  9  25,  and  I  Z  25  for  the  cases  of  Q  =50 
and  ICO,  with  =  1.25  and/^-^=  0.1  .  A  more  complete  coverage  of  possible 

values  of  ^  and/^-i  are  planned  in  the  near  future.  Realistic  values 
of  ^  cover  a  small  range  so  that  1,25  is  taken  as  representative.  These 
represent  the  physical  limits  of  the  parameters  and  any  meaningful  solution 
should  be  within  fhese  limiis. 

In  the  results  to  date  for  each  of  the  two  cases  discussed  above, 
the  line  of  intersection  of  the  F  =  0  and  the  G  =  0  surfaces  did  not 
intersect  with  the  H  =  0  surface.  With  the  higher  value  of  0  ,  the 


surfaces  are  further  from  simultaneous  solution.  As  yet,  tha  affect  of 
on  tfie  solution  has  not  been  observed.  The  possibility  that  a 
relaxation  or  modification  ot  the  Chapman-Jouget  condition  may  be  necessary 
in  order  to  obtain  a  simultaneous  solution  Is  presently  being  Investigated, 


APPTOIX  b 


System  of  Equations  end  Variables 


Variables 


f  =  PVp. 

f  ' 

2  = 


Equations  (!)  Hugoniot  relation  across  disturbance 

(2)  Isoenergetic  relation  between  disturbances 

(3)  Isentropic  relation  between  disturbances 
(6)  Chapman- Jouget  cond i tion 

These  ere  easily  reduced  to  three  equations  In  three  variables 

5^  ,  ^  ,  and  t2  ,  by  use  of  the  Isentropic  relation  to  ellrnlnate 

r .  Equations  (1),  (2)  and  (3)  are  replaced  by  Equations  (4)  and  (5J 

which  have  no  dependence, on  ^  .  The  terms  and y>j^  will  be  related 

to  ^  and  the  Injector  conditions  by  the  assumption  Vi-hlh.  and  by 

^  Z- 

the  isentropic  relations  applied  in  the  unburned  gas  region.  The  system  of 


equations  becomes: 


r(<p,Q-,i)=0 

G,(cp^&,z)^0 


~  s  )  =  ^1-0 


r>. 


Thssc  equations  can  be  written  out  in  terms  of  the  foj  lowing  quanitltiesj 


The  four  factors  in  Equation  (6-1)  are  the  following: 


(5-2) 


Using  the  solution  for  0  ,  9  ,  2  ,  we  find  the  following; 


U-Jfc/  -  SE-i -  h 

^O-i  0-(l  ^  L 
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Figure  I 
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Nonlinear  optical  rocket  motor, 
spark  source  and  camera 
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Figure  2 


Schematic  of  optica!”  test  stand 
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First  chamber  section  with  pulse  gun  in  tangential  firing  position 


Figure  4 


Injector  plate  with  Z  velocity  probes  in  place  and  tangential  oriented 


P=  Pressure  transducer 
V=  Velocity  probe 


Early  schlieren  photograph  of  shock 

pulse 


Figure  5 


Precision  diaphragm 


Figure  6 


Pulse  gun  cutaway 


JPR  945 


Detail  of  velocity  probe 
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Figure  9 


Schematic  diagram  of  nonlinear  stability  limit 

tests 
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Figure  10 


Ffrsf  nonlinear  combustion  instability  model 
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